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MORSAIR FIT DATA=GAMMA SILICON K=FACTOR, 
THERHONUCL SOURCE IN REAL ALR AT Bde KM. 
SAMP.ING ALTITUDES 80451005120, AND 143 KM. 
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SECTION I 
INTKODUCTION 


This report presents the results of fifteen computer 
runs made with the Kaman Sciences MORSAIR Monte Carlo computer 
program. These calculations were performed for the Air Force 
S Weapons Laboratory (AFWL/SAT) under the direction of Capt. 
| Raymond A. Shulstad. 


Ee? The objective of this effort was to provide AFWL with a 
set of "real" two-dimensional atmospheric neutron and seccndary 
z gamma transport data for use in assessing the adequacy of mass 
scaled uniform air caiculations. Previous studies at Kaman 
(ref. 1) and elsewhere (refs. 2, 3, 4) have indicated that 
large inaccuracies can result from scaling at high altitudes. 
The results presented in this report show the extent of these 
inaccuracies for two unclassified source spectra for both 


neutron and secondary gamma doses at al: _tudes from 5 to 80 
kilometers. 


For each MORSAIR run, the neutron and secondary gamma 
doses in silicon and tissue as well as the total particle 
fluence were calculated using a Monte Carlo technique at 
more than 150 detector locations about a point isotropic 
source in the atmosphere. The spectra used in these runs 
were an unclassified fission source and an unclassified 
thermonuclear source, both of which were provided by AFWL 
(ref. 5). In addition to a run in homogeneous air with 
each source spectrum for the purposes of verifying the 
tracking and scoring techniques of the code, eight runs with 
the fission source and five with the thermonuclear source 
were made in the variable density atmosphere. Section II of 


this report describes in detail the energy bins, spectra, 
altitudes and dose response functions used in the calculations. 
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In section III the homogeneous atmosphere neutron and 
secondary gamma transport data supplied by AFWL and used in 
this study are described. The two homogeneous atmosphere 


MORSAIR runs made for checkout purposes are discussed. 
Also, in section III, the statistical uncertainties of the 
t Monte Carlo data, the method used to fit these data, and the 
accuracy of these fits are briefly described. Finally, the 


method used to determine the K-factor, defined as the ratio 
2 


of the 4nR? dose in real variable density air to the 47R 
dose in homogeneous air, is described. 


In section IV, a brief discussion of some of the more 
obvious trends of the results is presented. The qualitative 


dependence of the K-factor on the source and detector altitude, 
| mass range, source spectrum, and dose response function is 

| described. The two primary effects which cause deviations 

| from the scaled homogeneous air results, namely a leakage 
effect and a mass distribution effect are also described. 


A summary of the calculations performed for this study 
and the conclusions which may be drawn from these results are 
presented in section V. 


In Appendix A, the tabulated fit coefficients of the 
one-dimensional ANISN 4nR? dose data (used in defining the 
K-factors) are shown for the two source spectra. All of the 
4nR° doses have been plotted for both neutrons and secondary 
gammas and these are included in Appendix A. In Appendix B, the 


coefficients for the fits to the MORSAIR two-dimensional data 
are shown for each sampling altitude, dose response function, 
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and particle type (neutron or secondary gamma). In Appendix C, 
plots of the fitted silicon dose data and the silicon K-factors 
for both neutrons and secondary gammas are shown for all 15 
MORSAIR runs. 


Be OC 18 


ee 
wg 
+ “Pe 
oe 
3 * 
vy 
r, 
cee 
aa 
a =f 
ae 
as 
we 
Soke 
ee, 
< aii 


SECTION II 
CALCULATIONS 


1. THE MORSAIR DRUGRAM 


The Kaman Sciences MORSAIR program is an extensively 
modified version of the ORNL MORSE program (ref. 6) developed 
specifically for Monte Carlo calculations of radiation 
environments in a variable density atmosphere based on the 
1962 Standard Atmosphere model (ref. 7). The multigroup 
cross section module of MORSE was used without modification, 
but the geometry and random walk modules were revised and a 
new scoring routine using an "extended path" or expectation 
boundary crossing estimator for scoring in concentric annular 
rings was added. In previous studies, the MORSAIR code has 
been used for calculating radiation environments at source 
altitudes from 20 to 65 kilometers. The results of those 
calculations have been checked against other Kaman real air 
transport codes, and codes from other agencies (ref. 8). It 
has been found to be a useful code for predicting radiation 
environments .t high altitudes where traditional methods of 
scaling uniform atmosphere results are inadequate. 


The latest version of MORSAIR was written for use on 
the CDC 7600 computer and utilizes the fast access large 
core capability of that machine so that very large problems 
in terms of storage requirements can be efficiently handled 
by the code. In this study the time integrated neutron and 
secondary gamma doses and associated standard deviations 
were Calculated for 58 energy groups at as many as 190 
detectors spaced about the source. The differential energy 
spectra were recorded on magnetic tape so that doses other 


than those shown in this report can be calculated. 
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Earth curvature effects were not included in the 


calculations, because these effects have been shown to be 


aw. 


negligible at the altitudes and ranges of the calculations 
performed for this study (ref. 9). 


2. INPUT PARAMETERS 


The air cross sections used in the MORSAIR runs were 
the DLC-31 multigroup set distributed by the Radiation Shielding 
Ei ¢ Information Center at ORNL (ref. 10). These cross sections 


| were prepared from ENDF/B-IV data using a fission spectrum 
weighting function. A third order (P3) Legendre expansion 
was used to represent the angular variation of the cross 
sections. The nitrogen and oxygen cross sections were mixed 
to form the macroscopic cross sections of air with the 
following composition (ref. 5): 


Density: Leli mg/em* 

Volume percentages: 79 percent nitrogen 
21 percent oxygen 

Number densities: 3.6609 x 1029 nitrogen atoms/cm* 
9.7316 x 1078 oxygen atoms/cm 


A copy of these cross sections was provided to Kaman by 

AFWL. The 37 neutron and 21 secondary gamma energy groups of 
the DLC-31 cross section set are shown in table 1. In table 

2, the fraction of source neutrons in each of the groups is 
shown for the two unclassified spectra used in the calculations. 


Three different response functions were used to weight 
the MORSAIR energy spectrum at each detector. Table 3 shows 
the neutron silicon and tissue dose response functions, and 
table 4 shows these response functions for the secondary 
gamma groups. The response functions shown in tables 3 and 
4 were provided by AFWL (ref. 5). In addition to the silicon 
and tissue doses, the total number fluence was also calculated 
for both neutrons and secondary gammas. 
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TABLE £ 


WoNeSsL ENSAGY GROUP STRUCTURE 


NEUTQON ENERGY (MEV) 


1.-9O4UE 40 1 
L.6905E¢02 
Le49LSE+I1 
1.46191F¢91 
1.-394UE4)1 
L294 E401 
Le2214E +91 
Le18S2E O01 
1.099E*]1 
220468454)0 
B.1A8738E499 
7 4082E +06 
66 3763E+039 
4. 9H5%E 400 
&e7237E409 
G0657EFIS 
3.0119E 400 
2.3452E+00 
2.3069F+09 
1.8208£4609 
1.1080E€¢00 
5.5023E-94 
1.5764E-94 
L1.1109E-d1 
5 62475E-92 
2e4788E=-)2 
2e1875E-02 
1.033%E+92 
Se 35hor-)3 
Le2341E 993 
5 eo BOI4E S04 
1.0139 E-04 
2.9023E-05 
1,.0677£-295 
3-0590E-05 
£e1254E 205 
4o1490E<97 


LeASUSECOL 
L.49L8E¢01 
1e419LE 601 
Le 3B40E 401 
1e2Q4K0E401 
1.2214E+01 
1.e1052E¢01 
1-0000E+01 
9.174546 400 
4218736400 
7.40826 ¢00 
&.3763E400 
4. 9659EF 400 
4.7237E%00 
&.0657F +00 
S.9LL9E 400 
223852E+00 
22 3059E 400 
16326396400 
1.1080E+00 
50 5023E201 
1¢5764E=61 
1.-1109E-01 
5 e2475E=02 
204 78B8E M02 
2.1375E-02 
1.0333t-02 
32354+6E-03 
1.23541E-03 
5 eo B294F 904% 
LeJLZOE-C4 
2.9023E-05 
1.0677E-05 
3eISIDE-06 
161254E-°06 
4.e1400E-07 
1.0000F-11 
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GAMMA ENERGY 


LeGQOOE+OL 
L.0009E+O4 
8.9000+00 
7.0900E+00 
6.0000E+00 
5 0000E +00 
4. 0000E+00 
3.2000€¢00 
2250008400 
2.0000€+00 
13000400 
1.0000E+00 
7.00906=-01 
4 5000E=01 
3.0000E=01 
1.5000E-01 
1.00006-91 
7.0000E-02 
4S Q00E=02 
3.0000E-02 
2.9000E-02 


(MEV) 


1.,0000E+01 
8.0000E+99 
7, 0000E+00 
6.0000E+00 
5.0000E+00 
4. 0000E+00 
3.0000E+00 
2.5000E+00 
2.0000E+00 
1.50006+00 
1.0000E+00 
7.000NE-01 
4. 5000E-01 
3.0000E-01 
1.5000E-01 
1.0000E-01 
7.0000€-02 
4e5000E-02 
3.0000€-02 
2.0000£-02 
1.00006-92 
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TABLE 2 


NEUTRON SOURCES 


NEUTRON ENERGY (MEV) 


1.9640£+01 
1.-6905E 404 
Le4dLBE+O1L 
LeGLILE COL 
L.3840E+01 
1.2354NE+O1 
1.22thE +34 
1.1052E +01 
LeO900E+01 
920484E+00 
BeLX73E+00 
7.4632E 4090 
6-3763E4+00 
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2e3069E+00 
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4 .9659E+00 
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4 .9657£ 400 
3.0419£+00 
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2.3059E+400 
1.4268E400 
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2e47IKBE“N2 
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323546E°03 
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5 eo d294E-04 
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TABLE 3 


NEUTRON OSE RESPONSE FUNCTIONS 


NEUTRON ENERGY (MEV) 


1.9640 E491 
1.6905E+91 
L.G4914E+91 
Le419te+I1 
1.3840E¢9 t 
1.284NE tot 
Le2214E+dL 
LeLO52E+IL 
L.O000E+I1 
9. 0484E+90 
AL LA73EFO0 
7.408200 
663763E+99 
4 .9659E+00 
ee 7237E+00 
4 0OSTE+I0 
SeOLLYE400 
2038522 +00 
2.3069£+90 
1682686420 
L.1080E+00 
5eS5U23E-91 
1.65754E*)1 
1.1109F -01 
5 e2475E-92 
2047ARE=22 
21 475E=92 
1.03336-02 
3.35465-03 
1.23416-93 
5 B294E-04 
1.0139E-04 
2.9023E=95 
1.0677£-05 
3.0590E-05 
1.1254E-05 
4e1400E=a7 
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1.9905E+01 
1.4918E+01 
Lo4L9LE+01 
1.3840E¢01 
1e2840E401 
1.221L4E401 
1.10526+01 
1.9000E+01 
2 e0494E 400 
3e1873E+00 
7.%032E+00 
6.37656 +00 
4 2 3659E+00 
4e7237E+09 
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2-3852E+00 
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1-8268E+00 
1.1030EF+00 
5.5023E-01 
165764E*01 
1-1109F-01 
5 e2475E-02 
2,.4738E=*02 
2e1875E=02 
LeUS33E-02 
323546F=-03 
1.2341E=03 
5 B294E=“04 
LeD1L30E-04 
209023E-05 
1.0677E=05 
3-9590E-06 
1.1254E-06 
Gel4DDE-N7 
1.9000b-11 


TISSUE OSE 
RAD/ (N/C42) 
8.6724E=+09 
7.4190E-09 
6e8115E-09 
§e5447E°09 
6e1473E-09 
5¢954BE-09 
5 e893EE=09 
5.5508E-09 
5 e2882E-09 
520473E-09 
5 .0045E-09 
%o7595E209 
4e483LE 209 
4e2531E-99 
#e1711E-03 
$e9784E209 
3.3905&-09 
3.13776-99 
3.0345E-09 
2.6393E-09 
200570E-09 
1.3330E-09 
7.6228r-10 
504890E"19 
3.1164E-10 
2.0739E-10 
1.4662E-10 
§e6143E"11 
2e275B8E-11 
90131512 
3.6632E-12 
1.1759E-12 
1.1095€-12 
1.6117 E-12 
2e7T415E-12 
445 70E-12 
Le1i238E-1i1 
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SILICON DOSE 


RAD/ CN/CM2) 
1.9106E-09 
1.7792E-09 
1.6813E-09 
1.6231£-09 
1.51-44E-09 
1.3851E-09 
1.237 0E-09 
1.0530€-09 
8.7897E~10 
7.9629E-10 
7.8141E-10 
&.7092E=10 
2.1394E-10 
1.826 7E-10 
104195E=10 
1.058 2E-10 
1.0006E-40 
6.2995E-11 
9.477 BE-14 
6.5328E-11 
4.9785E-11 
3.1515E-11 
1.7897E-12 
2.8022E-12 
1.232 7E-12 
7.908 4E-13 
5.893 06-13 
2.980 4E-13 
1.0498E-13 
4 3305E=14 
1e4421E=14 
4 SH95E=15 
3.937 76-45 
5.6286E-15 
90402 3E=15 
105390E=14 
7e4244E=13 
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TABLE 4% 


GAMMA DOSE RESPONSE FUNCTIONS 


GAMMA ENERGY 


1-4000E+9 1 
1.0000E+01 
800006409 
7.090CE+00 
6-0000E+00 
5.000)€ +400 
4.9000E+00 
3.0G05E+00 
2-50 0NE +90 
2.0C0fF +09 
1.500°E+00 
1.000nE+00 
7.-0000E-01 
&SIODE“91L 
3-O0690E-I1 
L1.500vE-01 
1.000j3E-91 
7.0003E-02 
4 S003 02 
3.009.5-02 
2U009E-02 


(Mev) 


1.0000E401 
8. J000E+00 
7.9000E+00 
660000400 
5. 0000£+00 
4. 0000F+00 
3.20006¢00 
2.500400 
2.0000E+00 
1.35000E+00 
4000006400 
7.2090E-014 
4 eSOQ0E-01 
3.00006-01 
1.65090E*01 
1.9090F-014 
7.9030E-02 
445000602 
$.9000E-02 
2.20006-02 
1.1000E-02 


TISSUE DOSE 
RAN (N/CM2) 
207431E=09 
2025E4E-09 
169840E-09 
1.7922E-09 
£.5928E-09 
1.3897E-19 
1.1803E-09 
1.9098E-39 
3.5320¢-10 
7.4281E-10 
5.30 306-10 
4e2393E-10 
269695E-10 
1.9283E-10 
100770E-10 
0 9383E=11 
304 315E=11 
2.9479E-11 
4037506-11 
306647E-11 
3.2504E-10 


SILICON DOSE 


RAD/(N/CM2) 
32418 GE-09 
2e5712E-C69 
2.1612E-09 
1.8991€-09 
16636 7E-09 
123835E-09 
101335E-09 
Go4SS4E-10 
8203 4E-10 
&e5343E-10 
5 e2B846E-10 
3.8505E-10 
2e7122E-10 
1.777 2E-10 
12045 9E-10 
7.985 4E-11 
1.454 3-10 
3e44141E 210 
8o2679E-10 
20649 3E-09 
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3. DESCRIPTION OF RUNS 


A summary of the 15 MORSAIR runs made for this study is 
shown in table 5. In addition to the homogeneous atmosphere 
runs, eight runs were made using the fission source in real 
air at altitudes from 5 to 80 kilometers, and five thermo- 
nuclear source runs were made at altitudes from 20 to 80 
kilometers. A minimum of 50,000 initial neutron histories 
were followed for each run, and the secondary gamma production 
rate was adjusted so that approximately the same number of 
gammas was produced in each run. For the lower altitude 
and homogeneous runs, the computer time required for 50,000 
histories was greater than 30 minutes of CDC 7600 execution 
time. At the higher altitudes an appreciable number of the 
neutrons and secondary gammas escape out the "top" of the 
atmosphere after a few scatterings so that as many as 150,000 
initial neutrons and approximately the same number of sec- 
ondary gammas could be followed in less than 30 minutes. 


A time cutoff of 20 seconds was used for each neutron 
history. This ensured that almost no neutrons were lost 
because of this cutoff. 


The atmosphere model used in MORSAIR provides for a 
continuously varying density from sea level to 200 kilometers 
and very closely approximates the 1962 Standard Atmosphere 
model (ref. 7). Above 200 kilometers, a void is assumed and 
particles reaching this altitude are allowed to escape. The 
lateral extent of the atmosphere varied with each run to 
provide a large buffer beyond the last detectors of interest 
at each sampling altitude. This ensured that lateral leakage 
effects on the results would be negligible. Any particle 


reaching the ground (which only a very few did) was terminated. 
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Particle histories were terminated primarily because they 
escaped from the top of the atmosphere or their weight 
became so small that they could not contribute a significant 
| score at any detector. 


SECTION III 
RESULTS 


l. HOMOGENEOUS ONE-DIMENSIONAL TRANSPORT DATA 


The one-dimensional homogeneous air neutron and secondary 
3 gamma transport data used as a basis of comparison in this 


report were provided by AFWL (ref. 12). These data were 


¥ ; generated by appropriate source spectrum weighting of eat 
a i fits to Burgio's ANISN results (ref. 11). The resulting 47R 
* dose data were then fit to the equation: i 


ln [4nr? Dose | = A+ Bx + Cx? + px2/2 + Ext /2 1/3 


+ Fx + Ginx 


for 0.1 < x < 300 gm/cm’ 
( 
BT i where x is the mass range‘ of air or areal density in gm/om?. 
: The seven coefficients representing the fits to the neutron ' 
and secondary gamma homogeneous atmosphere data are tabulated 
in Appendix A for the two spectra and three dose response - 
e functions used in this study. 
i For very small areal densities, between 0.0 and 0.1 
: 1 gm/em? of air which are of interest particularly at very 
Ss ’ high le eke a linear SRE on ace oP between the a eee 
eae ANISN 47R” dose at 0.1 gm/cm' and the uncollided 47R” dose 
rat (i.e., the sum over all energy groups of the products of the 
Bat source fractions and dose response functions) at zero mass 
2 ¢ range was used to find the 4nR? neutron doses. This method 
& ¥ is commonly used in mass integral scaling codes to define 
Mex atmospheric neutron environments at small mass ranges. A 
- similar interpolation technique tar used for the secondary ‘ 
he gammas between zero and 0.1 gm/cm™ of air, except in this : 
o3 case a 4nR? dose value of zero was used at zero mass range. 7 
an *See footnote on page 35 for a definition of mass range. 
*™ 
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The same energy groups, cross sections, sources, and 
dose response functions described earlier for the MORSAIR 
runs were used in the uniform air calculations. 


| : 2. MORSAIR VERIFICATION IN HOMOGENEOUS AIR 


® An option exists in the MORSAIR code to replace the 
variable density atmosphere model with a homogeneous atmosphere 


Be i model. As an initial check on the scoring and tracking 
techniques employed in the code, one run for both the thermo- 
nuclear and the fission source spectrum was made using this 
* option. It was anticipated that the results of such a run 
would agree with the one-dimensional ANISN infinite air 
results provided by AFWL from Burgio's calculations using 
the same spectra and response functions. Figures C~-1 and 
C-3 of Appendix C show that for all sampling altitudes the 


. difference between the fits of the ANISN homogeneous data 
and fits of the two-dimensional MORSAIR homogeneous data is 
. less than 15 percent out to 100 gm/om? and within 25 percent 


out to 200 gm/om? for a fission source. Similar results for 
a thermonuclear source are shown in figures C-49 and C-51 of 
Appendix C. It should be emphasized that the solid lines 
shown in these and other figures of Appendix C represent the 
fitted one-dimensional homogeneous air ANISN d. a4, and the points 
are calculated from the fits of the MORSAIR data. It should 
also be pointed out that in the figures immediately following 
those mentioned above, namely figures C-2, C-4, C-50, and 
C~52 of Appendix C, the K-factors presented are actually ratios of 
fits of the homogeneous MORSAIR data to the ANISN data. All 
of the data in Appendix C is presented in this manner, i.e., the 
fitted homogeneous and real air data are shown first followed 
by the K-factors calculated from this data for each MORSAIR 
< run. The differences between the ANISN and MORSAIR results 

for a homogeneous atmosphere are within the statistical 
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uncertainties of the MORSAIR data for these runs, and it was 
concluded that the MORSAIR tracking and scoring methods were 
correct. 


3. STATISTICAL UNCERTAINTIES IN THE MORSAIR DATA ” 


At the lowest source altitudes, the 50,000 neutron and 
secondary gamma histories resulted in standard deviations in 
the total doses of less than 15 percent out to 100 gm/em*, and 


less than 25 percent to 200 gm/em” in most cases. At source 
altitudes above 20 kilometers, the expectation boundary 
crossing scoring method is more efficient, and the statis- 
z tical uncertainties in the data were somewhat less. 


4. FITS TO THE MORSAIR DATA 


To smooth the MORSAIR Monte Carlo results and to reduce 
the quantity of data to a more usable form consistent with 
the one-dimensional results, all of the 4m? doses were fit 
to the same seven parameter equation in areal density 
described above. For the MORSAIR data, it was found that 
dropping the last term, G ln x, did not affect the accuracy 
of the fits significantly and in some cases resulted ina 
slight improvement in the fit. For this reason, only six 
coefficients are used in the fits to the MORSAIR data shown 
in Appendix B. 


The technique used to fit the data was a standard linear 
weighted least squares technique such as described in reference 
13. The fitting method, it should be noted, was applied to the 
logarithm of the 4nR? dose data where each data point was given 
a weight inversely proportional to the square of the logarithmic 
uncertainty in the MORSAIR data point. Any data point with a 
standard deviation greater than 50 percent was deleted before 


Sy : the fit was applied. The number of points removed for this 


3 
an ES reason was small for areal densities less than 200 gm/om. 
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The fit coefficients for each run are tabulated in Appendix 
B. tn addition, the range (in terms of minimum and maximum 
areal densities) over which the fit applies ig shown. The last 
column in the tables, "RMS PCT DIFF", is the weighted root mean 
square percentage difference between the fit data and the 
actual MORSAIR data, i.e., 


"RMS PCT DIFF" 


AnR? actual MORSAIR dose 


2 


where Y, = In 
4nR” £it MORSAIR dose 


and the weight, Wye was defined as: 
-2 


W. 
1 


ln (1 + 3) 


where Os fractional standard deviation of the 
MORSAIR dose. 


As an example of the adequacy of the six parameter fit 
to represent the data, figure 1 shows both the actual MORSAIR 
data points and the values calculated from the fit for 
coaltitude neutron silicon doses from a thermonuclear source 
at 20 and 40 kilometers. Figure Z shows the MORSAIR data 
points and the fit for the secondary gamma silicon dose. In 
general, it was found that this fitting technique resulted 
in a reasonably good representation of the actual data. 


Ds THE K-FACTOR 


The adequacy of mass integral scaling of uniform air 


results to define real air environments can be conveniently 
described in terms of a K-factor which was first used by 
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Figure 1. MORSAIR Data Points and the Fit to the Data Points 
for 41R? Neutron Silicon Dose, Coaltitude Sampling 
from 20 and 40 km Thermonuclear Sources 
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Marcum in his early high altitude transport studies (refs. 
2, 3). The K-factor as used in this report is defined as 


The 4mR? dose at a detector as calculated in a 
"real" variable density atmosphere (from fits to 
MORSAIR data in this case) 


The 4nR° dose received at the detector as cal- 
culated in homogeneous one-dimensional infinite 
air (from fits to the ANISN data in this case) 


K-factor = 


The K~factor, then, is a direct measure of the error 
associated with mass integral scaling of infinite air data 
due to the variable density nature of the atmosphere and can 
be used as multiplicative correction factor to the scaled 
data. A K~factor of 0.5, for example, indicates that in the 
real atmosphere the dose received is only one half of the 
dose that would be predicted by scaling methods. 


Marcum's early studies and more recently those of Keith 
(ref. 1) and cthers (refs. 14, 15) have shown that the neutron 
K-factor is a function of source altitude, detector altitude, . 
slant range, source energy, and the dose respcnse function. 

The neutron K-factor, it has been shown, can vary from less 
than .1 at high source altitudes to greater than 5 for 
detectors below the source and at large areal densities. It 
is believed that while this study is more extensive in terms 


of the number of calculations performed and the range of 
altitudes and mass ranges covered, the data trends shown in 
this report for neutrons are in substantial agreement with 
the earlier data of Marcum and Keith. No previous attempts 
to calculate the transport of secondary gammas in a real 


atmosphere have been published. 


Although there are many situations of practical interest 
in which the K-factor is near unity so that scaling results 
in small errors, there are also situations in which large 
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It is hoped that the data presented in this report will 
4 ; allow a user to recognize those situations where scaling may 
be inaccurate. 


4 | errors can result from neglecting variable density effects. 


a | ; In Appendix C, the AmR? silicon dose and the silicon K- 
‘% j factor are plotted as a function of areal density for each 
of the MORSAIR runs. The trends of the silicon dose K-factors, 
yo as will be shown in section IV, are similar to those of the 
| tissue dose and the total fluence K-factors so that in the 


interest of brevity, only the silicon dose and K~factor plots 
| are shown. 


. 


Definition of Mass Range (RHOR): 


In homogeneous air, 


RHOR = p(zs)R 


‘and in 2-D real air, 


is o(z)dR 
0 


density as a function of altitude; 


RHOR 


where 


p (2) 


en 


p (zs) = density at source altitude (zs); 


oh: 


ae. 


R = slant range between the source and a receiver point. 
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SECTION IV 
DISCUSSION 


The K-factor, as defined previously, can be considered 
as a multiplicative correction factor to be applied to the 
scaled homogeneous air dose to account for the transport of 
neutrons and secondary gammas in the real variable density 
atmosphere. Since many computer codes used to predict neutron 
and secondary gamma environments in the atmosphere (Such as 
ATR and SMAUG) use the method of mass integral scaling of 
homogeneous air transport data, it is believed that presenting 
the MORSAIR results in terms of K-factors allows a direct 
means of assessing the accuracy of these codes for source 
spectra similar to those used in this study. 


The results of this study show that the K-factor is a 
complicated function of several parameters. In this section 
the dependence of the K-factor and therefore the validity of 
mass integral scaling of homogeneous air data on these 7 
parameters is described in a qualitative manner. A detailed 
explanation of the differences in the transport of neutrons 
and secondary gammas in homogeneous air and in a variable 
density atmosphere was beyond the scope of this effort. 


l. K-FACTOR DEPENDENCE ON ALTITUDE AND MASS RANGE 


In general as the source or detector altitude increases, 
more particles escape out the top of the atmosphere. Therefore, 
when compared to the mass integral scaled dose, the dose in 
a real atmosphere can be substantially less. This leakage 
effect has been demonstrated in many high altitude transport " 
calculations. Pigure 3 shows the coaltitude neutron K-factors 
for silicon at several source altitudes as calculated by MORSAIR . 
for a fission source. For a detector at 40 gm/om? from the 


source, the neutron K-factor is nearly unity for sources below 
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Figure 3. Neutron Coaltitude Silicon K~-factors at Several 
Source Altitudes for a Fission Source 


20 kilometers, decreases to about 0.5 for a source at 30 
kilometers, and is less than 0.2 for a 40 kilometer source 


cece nme etn tr alba -thece * 
a aa rameters 


altitude. The secondary gamma silicon K-factors exhibit a 
similar behavior as shown in figure 4, except the K-factors 
fall off even more at the higher source altitudes. Similar 


trends in the K-factors occur if the source altitude remains 
constant and the detector altitude changes as shown in 
figures C-27 and C~29 of Appendix C, for example. 


In addition to this expected decrease with source and 
detector altitude due to leakage effects, it is also obvious 
from figures 3 and 4 that the K-factors are a strong function 


of the amount of air (areal density) between the source and 
detector. For detectors coaltitude with the source, the K- 
factor at high altitudes decreases monotonically as the 
areal density increases, but at the 20 and .0 kilometer 
source altitudes, the K-factor curve is more complicated. 


Figure 3 shows that for these source altitudes, the 
neutron K-factor initially decreases with mass range as it 
does at higher altitudes due to leakage of particles from 
the atmosphere. At mass ranges of 80 to 100 gm/cm?, however, 
the neutron K-factor begins to increase and at very large mass 
ranges it can be greater than 4 indicating that at these 
large ranges the actual dose can be more than four times the 
dose that is predicted by scaling the homogeneous transport 
data. This dose enhancement has been observed in earlier 
results (ref. 1), and is apparently due to the fact that a 
substantial number of the neutrons which arrive at these 
detectors at large mass ranges travel large distances in the 


lower density air above the detector and thus traverse less = 
air than is predicted by the simple straight line mass integral 


between the source and the detector. They are therefore 
attenuated less than what would be predicted by the scaled 
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Figure 4, Secondary Gamma Coaltitude Silicon K-factors at 
Several Source Altitudes for a Fission Source 
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homogeneous atmosphere results. This effect has been 
referred to as a “short circuiting" (ref. 15) or "mass 
@istribution" (ref. 14) effect. While it is most apparent 
for detectors coaltitude and below the source at large mass 
ranges, this effect also occurs to a lesser extent for 
detectors above the source at large distances as is apparent 
from figure C-18, for example. 


2. K-FACTOR DEPENDENCE ON DOSE RESPONSE FUNCTIONS 


The reduction or enhancement of the neutron dose in the 
real variable density atmosphere over the mass integral scaled 
dose depends not only on the source and detector altitudes 
as shown previously, but also on the dose response function 
used. For a source at 20 kilometers, it is shown in figure 
5 that the neutron K-factors for the silicon and tissue 
doses and the total fluence are similar for sampling altitudes 
coaltitude with the source. At 40 kilometers, however, figure 
5 shows that the neutron K-factors can be quite different, 
with the tissue dose K-factor being only half as large as the 
silicon K-factor, and the K-factor for fluence even smaller. 
These differences in the K-factors at high altitude are 
related to the fact that the response functions are different 
and the energy spectrum arriving at a detector in the real 
atmosphere can be quite different from the spectrum calculated 
in homogeneous air at the same areal density. From this 
example it is clear that the neutron real air effects depend 
strongly on the dose response function used to weight the 
neutron fluence. 


While the secondary gamma K-factors are strong functions 
of source altitude and mass range from the source as shown in 
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figure 6, they are not so dependent on the response functions 


as was the case for neutrons. 


ce K-FACTOR DEPENDENCE ON SOURCE SPECTRUM 


In figure 7 the coaltitude neutron silicon dose K-factor 
is shown for both the thermonuclear and fission source 
spectra at source altitudes of 20 and 40 kilometers. At 20 
kilometers, the differences between the two K-factors are 
small at all areal densities, despite the fact that the two 
initial emitted spectra are quite different. At 40 kilometers, 
however, the Silicon dose K-factors for the thermonuclear 
source are much larger than the K-factors for the fission 
source at areal densities less than 70 gm/em-. At 40 gm/em" , 
for example, the coaltitude neutron silicon K-factor for a 
thermonuclear source spectrum is twice as large as the K- 
factor for a fission source. Similar results for the 
secondary gamma doses are shown in figure 8. At 60 gm/om? 
the silicon K-factor for secondary gammas from a thermonuclear 
source is less than one-third of the silicon dose K-factor 
from a fission source spectrum coaltitude from a source at 
40 kilometers. At 20 kilometers, however, the secondary 
gamma K-factors are nearly identical. 


It is apparent from figures 7 and 8 that mass integral 
scaling can lead to very large errors in the neutron and 
secondary gamma environments in the atmosphere, and that the 
magnitude of these errors depends rather strongly on the 
source spectrum. In general, the K-factors for the fission 
source are less than those for a thermonuclear source indicating 
the scaling errors for this source spectrum are even greater 
than they are for a thermonuclear source. 
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Figure 6. Silicon, Tissue, and Fluence Secondary Gamma K-factors 
Coaltitude from a 20 and 40 km Thermonuclear Source 
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Figure 7. Coaltitude Neutron Silicon Dose K-factors at 20 and 40 
km from both a Fission and Thermonuclear Source 
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Figure 8. Coaltitude Secondary Gamma Silicon Dose K-factors at 
20 and 40 km from both a Fission and Thermonuclear 
Source 


4. SECONDARY GAMMA K-FACTORS AT SMALL MASS RANGES 


There is one feature of the results for secondary gamma 

doses at the highest altitudes, 60 kilometers and above, 

that deserves an additional comment. Figure C-84 in Appendix 

C is a good example of this data for the case of a thermo- 

nuclear source at 80 kilometers. The silicon K-factors for 

the higher sampling altitudes are very large at low mass 

ranges. The K-factor, as was pointed out in section 3, is 

Bi the ratio of the 4mR? dose in real air to the 4mR* dose in 

homogeneous air. The homogeneous results for this study 

ao were based on ANISN code calculations. Because of the 
uncertainties in these results at small mass ranges, most 
scaling codes use an interpolation method between 0.0 and 

about 1.0 gm/om?. Such a method was used to produce the 

: homogeneous data shown in Appendix A. No such approximations , 

- or assumptions were necessary in the MORSAIR results. So it 

a should be noted that while the K-factors at these small mass 

ranges may not necessarily reflect differences between the 

homogeneous and real air results, they do reflect differences 

between the real air results and those from typical scaling 


aS codes. 
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SECTION V 
. SUMMARY AND CONCLUSIONS 


In this report we have presented the results of 15 
‘MORSAIR Monte Carlo code calculations of neutron and secondary 
gamma transport in a real variable density atmosphere. The 
two source spectra used in these calculations were an unclassified 
fission spectrum and an unclassified thermonuclear spectrum. 
The source altitudes varied from 5 to 80 kilometers and the 
A energy dependent particle fluences were scored at more than 
150 detectors located about the source for each run. The 
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4 spectrum at each detector was weighted with fluence, silicon, and 
! tissue dose response functions. The 4nR? dose data has been 

‘- fit to a seven parameter equation in areal density at each 
sampling altitude of the MORSAIR run. Coefficients of these 

< fits and plots of the silicon dose and K-factors have been 
presented. We have shown how the K-factor, the correction 

to the scaled homogeneous transport data to account for 

variable density atmosphere effects, depends on several 
parameters. 


As stated at the outset, the objective of this study 
was to provide data which would lend some insight into the 
accuracy of mass integral scaling of homogeneous air transport 
data at highaltitudes. The data presented in this report 
demonstrate that mass integral scaling can lead to large 
errors particularly for source altitudes greater than 20 
kilometers. We believe these data can serve as the basis of 
' an algorithm for predicting the correction to the mass integral 
{ . scaled dose for sources and response functions similar to those 
used in this study. 
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APPENDIX A 


ANISN HOMOGENEOUS AIR DATA 


This appendix contains the coefficients of the fits to 
the neutron and secondary gamma transport ANISN code results 
provided by AFWL (ref. 12). The data was fit to the equation 


ln (40R? Dose) = A + Bx + Cx? + px3/2 


; : 2 
where x is the mass range in gm/cm”. 


1/2 1/3 


+ Ex + FX 


The coefficients of this fit for the different doses and 
sources are shown in Table 6. 


The 4nR? doses and fluences ascalculated using these 
fits for both neutrons and secondary gammas have been plotted 
t 

for the two source spectra in figures A-1 through A-6. 
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Figure A-1. Homogeneous ANISN Fit Data. 47R? 
Dose for a Fission and Thermonucle 


Neutron Silicon 
ar Source. 
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Homogeneous ANISN Fit Data. 47R? Neutron Tissue 
Dose for a Fission and Thermonuclear Source. 
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Figure A-3. Homogeneous ANISN Fit Data. 47R* Neutron Fluence 
for a Fission and Thermonuclear Source. 
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Figure A-4. Homogeneous ANISN Fit Data. 47R? Secondary Gamma 
Silicon Dose for a Fission and Thermonuclear Source. 
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Figure A-5. Homogeneous ANISN Fit Data. 4nR? Secondary Gamma 
Tissue Dose for a Fission and Thermonuclear Source. 
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APPENDIX B 
FITS OF MORSAIR REAL AIR DATA 


This appendix contains the results of a weighted least 
squares fit to the MORSAIR data. For each of the 15 MORSAIR 


runs, the Zollowing data is given at each sampling altitude: 
Ls The sampling altitude in kilometers. 
2. The dose response function. 


3. The minimum and maximum mass ranges at which the 
fit is considered valid. 


4. The coefficients to the fit equation: 


in (40R2 Dose) = A + Bx + Cx? + px3/2 + Exl/2 + Pxl/3 + Ginx 


where x is the mass range in gi/om? 


5% The root mean square percentage difference between 


the fit values and the actual MORSAIR data. 


These data are shown for both the neutron and secondary 
gammas. 
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APPENDIX C 
MORSAIR SILICON DOSES AND K-FACTORS 


For each of the 15 MORSAIR runs, this appendix contains 


the computer plotted results of: 


1. The fit of the MORSAIR ar R2 silicon dose versus 
the mass range ("RHOR" in the figures). The 
different point symbols were calculated from the 
fits to the MORSAIR data, and the solid line is 
the fit of the one-dimensional ANISN data. 


2. The silicon dose K-factors as a function of 


mass range at each sampling altitude. 


In order to reduce the number of plots, several sampling 


altitudes have been combined on a single figure. 
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FIGURZ Ceey 4ORSATA FIT GATA@4PIR*¥2 NEUTRON STI DOSE, 
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FIGURE C-50 YORSAIR FIT GATA*NEUTRUN SILICON K-FACTOR 
FISSION SOURC!, IN Rc AL AIR AT 0460 KM, 
SAMPLING ALTITUD=5 26 AND 4S KM. 
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FIGURE U-3¢2 AMURSAIR FIT DaTA-GaMna SILIGON K-FACTOR, 
FISSION SOURCE IN REAL AIR AT 60.0 <M. 
SAMPLING ALTITUDES 20 AND 4 KM. 
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FIGUR: 3°35 WMORSATR FIT DATA-4PIRF*2 GAMMA SI OSE, 
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FIGURE C-36 MORSAIR FIT OATA*GAMMA SILICON K=FACTOR, 
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FIGURE C*#37 MORSAIR FIT DATA~GPIR**2 NEUTRON SI DOSE, 
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FIGURE C-38 YORSAIRK FIT OATA-NEUTRON SILICON K-FACTOR 
FISSION SOURCE IN REAL AIR AT 60.0 KM. 
SAMPLING ALTITUDES 20 KILOMETERS. 
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FIGURE Ce$9 MGORSALF FIT DATA-4PIR**2 GAMMA SI UOSE, 
FISSLON SOURCE IN REAL AIR AT 84.3 KM. 
SAMPLING ALTITUOrS 20 KLLOMETCERS. 
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FIGURE 3°40 MORSAIR FIT UATA=GAMMA SILICON K*F ACTOR, 
FISSION SOGURGe IN REAL AIR AT 6060 KMe 
SAMPLING ALTITUDES 20 KILOMETERS. 
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FIGURe Ge4¥l YGRSAIR FIT OATA-4PLR**2 NEUTRUN SIL DOSE, 
FISSIGH SOURCK LN RIAL AIR AT 80.0 <M. 
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RUN NOo 9 FISSION ‘SOURCE 
SOURCE ALTITUDE= 80.000 KM NEUTRONS RADS-SILICON 


1 the 


m=40.0 KM 
QO KM 


YaPImRaw2 RADS-SILICON 


19°! 


0. 000 5.000 10, 000 20, 000 25 500. 30, 000 


15.006 
RHOR G/CM2 


FIGURE Ceui YORSAIR FIT OATA-4PIR**2 NEUTRON SI OSE, 
FISSION SOURCe LN RIAL AIR AT 80.0 <M. 
SAMPLING ALTITUDcS 40 ANO 60 KM. 


116 


~ 


re re TE TE er I Se 


ET 


RUN NO. 9 FISSION SOURCE 
SOURCE SLTITUDE= 80.000 KM NEUTRONS AADS-SILICON 


10 * 


m=40.0 KM 
@=60,0 KM 


K FACTOR RADS-SILICON 


10-4 
: 0.000 5.000 10,000 20.000 2000 30,000. 


15, 000 
RHOR G/CM2 


FIGURE C-42 MORSAIX FIT QDATA*NEUTRON SILICON K-FACTOR 
FISSION SOURCE IN REAL AIR AT 80.0 <Me 
SAMPLING ALTITUGES 40 AND 60 KM. 
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FIGURZ U-4S MORSale FIT UAIA*4PIR**2 GAMMA ST GOSE, 
FISSION SOURCE IN RLAL AIR AT 850d <M. 
SAMPLING ALTITUJES 46 ANU 6u KM. 
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FIGURE U-44 MOKSAIR FIT DATA-GAMHA SILICON K“F ACTOR, 


FISSION SOURCE IN REAL AIR AT 80.0 KM. 
SAMPLING ALTITUD<S 40 AND 69 KMe 
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FIGURE L=-46 4ORSAIR FIT DATA*NEUTRON SILICON K-FACTOR 
FISSION SOURCE IN REAL AIR AT 82.0 KMe 
SAMPLING ALTITUD:S 30410051204 AND 14d <M. 
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RUN NO. 10 THERMONUCLEAR SOURCE 
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FiGUR: Cebu MORSAIN FIT DATA-NEUTRON SILICON K-FACTOR 
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FIGURE C-f1 HORSALS FIT DATA-G4PIR*¥*2 GAMMA SI DOSE, 
THERMGNUCL SUURCE IN HOMO AIR AT 4.9 KMe 
ALL SAMPLING ALTITUDES 
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FIGURC C-52 MORSAIR FIT DATA-GAMMA SILICON K-FACTOR, 
THERMONUCL SOURCE 1N HOMO AIR AT 4.9 KM, 
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FIGURE Ce55 MURSAITR FLT DATA-GPIR*¥*2 GAMMA SI DOSE, 
THERMONUCL SOURCE IN RCAL AIR AT 20.0 KM. 
ALL SAMPLING ALTITUOCS 
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FIGUR~ C-56 MORSAIR FIT DATA-GAMMA SILICON K-FACTOR, 
THERMONUCL SOURCE IN RIZAL AIR AT 20.0 KM. 
ALL sAMPLING ALTITUDES 
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FIGURe 3°57 AORSAIR FIT OATA-GPIR¥*2 NEUTRON SI COSE, 
THERMOWUCL SOURCE IN REAL AIR AT 36-0 XM. 
ALL SAMPLING ALTITUDES 
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FIGURE C-v8 40KSAIR FIT DATA-NCUTRON SILICON K-FACTO? 
THERMONUCL SOURCE IN REAL AIR AT 30.0 KM. 
ALL SANPLING ALTITUD.'S 
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RUN NO. 12 THERMONUCLEAR SOURCE 
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Fe" GNUGL SUURCE IN RAL AIK AT 30.0 <M. 
ail SAMPLING ALTITUDES 


134 


meet 


¢ 
F 
I 


K FACTOR RADS-SILICON 


Ree 


i a Faey TCE RTE FORE — OE 
i | 


RUN NG. 12 THERMONUCLEAR SOURCE 
SOURCE ALTITUDE= 30.000 KM SEC GAMMAS RADS-SILICON 


10 ° 


m=i5.0 

@©=20.0 

¥ 4230.0 
a +240.0 
% =50,0 
46.0 


ZE8335 


0.000 4,000 6000 12000 16000 20.000 24,000. 
RHOR G/CM2 (X10) 


FIGU  ‘-uG -MORSALR FIT UATA-GAMMA SILICON K-F ACTOR, 
THERMONUCL SOURCE IM REAL AIR AT 300d KM. 
aL& SAMPLING ALTITUDZS 


135 


Ta eee FE Ey a Pe rs S re Fi 


Sere 


RUN NO. 13 THEAMONUCLEAR SOURCE 
SOURCE ALTITUDE= 40.000 KM NEUTRONS AADS-SILICON 
10" 
m=20.0 KM 
©=30.0 KM 
4240.0 KM 
+=60.0 KM 
=80.0 KM 


7" @=100.0 KM 
10 


10°? 


UxPImRae2 RADS-SILICON 


10°% 


i¢°" 
0.000 4. 000 8.000 12000 18000 20.000 24000 
RHOR G/CM2 (X10)! : 


FIGURc G61 YORSALR FIT UaTA-4PIR**2 NEUTRON SI DOSE, 
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FIGURe Cev6 MOKRSALIR FIT DATA=GAMMA SILICON K-F ACTOR, 
THERHONUCL SOURCE IN REAL AIR AT 69.0 KM. 
SAMPLING ALTITUD-=S 20 AND 4G KM. 
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FIGURE C72 MORSAIR FIT OATA=GAMMA SILICON K-FACTOR, 
THERMONUCL SOURCE IN REAL AIR AT 60.0 KM. 
SAMPLING ALTITUDES 605805100, AND 120 KM. 
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FIGURI 5°73 MOKSALR FIT OATA-4PI2**2 NEUTRON SI DOSE, 
THERMONUCL SOURCE IN REAL AIR AT 80-0 KM. 
SAMPLING ALTITUDES 20 KILOMETERS. 
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FIGURE C-75 MORSAIR FIT DATA-GPIR**2 GAMMA SI DOSE, 
THERMONUCL SOURCE IN REAL AIR AT 6U.G KM. 
SAMPLING ALTITUDES 29 KILOMETRE. 
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FIGURE C-76 YORSAIR FIT GATA-GAMMA SILICON K-F ACTOR, 
THERMONUCL SOURCE IN REAL AIR AT 8060 KM. 
SAMPLING ALTITUOD“S 20 KILOMETERS. 
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FIGURE G/B MORSAIR FIT OATA-NEUTRON SILICON K-FACTOR 
THERMONUCL SOURC:. IN RceAL AIR AT 6060 <Me 
SAMPLINC ALTITUDES 40 AND 60 KM. 


153 


bad ON sil a at a eet a aa aE a i a tla ik i ala Ra, em i ka ail Ml tea i eal maaan le OA ea aaa al, tala ih adil ed 
S 
? 


RUN NO. 15 THERMONUCLEAR - SOURCE 
SOURCE ALTITUDE= 60.000 KM SEC GAMMAS RADS-SILICON 


10°° : 


m=*40.0 KM 


®=60.0 KM 
| 
. 10 
: & 10 
oO 
Den 
' a | 
| Lo | 
i a 
| Ww 
Q 
ram 
o 
i N 4 
m 
” 
o 
m7 
—~ e 
a 
rm 
0" 
1 or? 
0.000 5.000 10.000 15.000 20. 000 25. 000 30,000 
RHOR G/CM2 is 


FIGURE C-/9 MORSALK FIT UATA-4PIR**2 GAMMA SI OOSE, 
THERMONUCL SOURCE IN RtAL AIR AT 80.0 KM. 
SAMPLING ALTITUDES 49 AND 60 KM. 
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FIGURE C-6d MOKSAIR FIT DATA-GAMAA SILICON K-FACTOR, 
TRERMONUCL SOURCE IN REAL AIR AT 80-0 KMe 
SAMPLING ALTITUDES 40 AND 60 KM. 
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FIGURE Codi MOKSAIK FIT DATA-4PIR*¥*2 NEUTRON SI OOSE, 
THERMONUCL SOURCE IN REAL AIR AT 80.0 KM. 
SAMPLING ALTITUDES 80910001205 AND 140 XM. 
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FIGURE C-862 MORSAIR FIT DATA=NEUTRON SILICON K-FACTOR 
THERMONUCL SOURCE IN REAL AIR AT 80.0 KM. 
SAMPLING ALTITUDsS 80410041205 AND 140 KM. 
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FIGURE C-83 MORSALR FIT DATA-4PIR**2 GAMMA SI DOSE, 
THERMONUCL SOURCE IN REAL AIR AT 806.0 KM. 
SAMPLING ALTITUOc%S 6041009120, AND 1460 KM. 
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FIGURE C-84 MORSAIR FIT DATA-GAMMA SILICON K-FACTOR, 
THERMOKUCL SOURCE IN REAL AIR AT 80.0 KMe 
SAMPLING ALTITUDES 8051005120, AND 149 KM. 


